Protected Rabi oscillation induced by natural interactions among physical qubits 
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For a system composed of nine qubits, we show that natural interactions among the qubits induce 
the time evolution that can be regarded, at discrete times, as the Rabi oscillation of a logical qubit. 
Neither fine tuning of the parameters nor switching of the interactions is necessary. Although 
straightforward application of quantum error correction fails, we propose a protocol by which the 
logical Rabi oscillation is protected against all single-qubit errors. The present method thus opens 
a simple and realistic way of protecting the unitary time evolution against noise. 
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Ch ! I. INTRODUCTION 

Decoherence of quantum states has been attracting much attention for long years Many methods have been 
proposed for defeating the decoherence. As compared with other methods [1, H, H], [H , quantum error correction (QEC) 
[1, 0, H, [HI has a great advantage of protecting against arbitrary errors if they only affect a single qubit (two-level 
system) in each logical qubit Q. Although QEC has been developed in the context of quantum computation, it is 
interesting and useful to apply QEC to protection of the unitary time evolution (Hamiltonian evolution) against noise. 

When trying to realize this, however, one encounters many physical problems, which are usually disregarded in 
discussions on the computational complexity [§}. For example, some physical process may be much more difficult 
to realize than another, even if the number of the necessary steps for them differs 'only by polynomial steps' [![. 
Furthermore, fabrication of a controlled-NOT gate, which is one of the elementary quantum gates, is very difficult 
because it requires fine tuning of the coupling constants of the interactions and high-precision switching of them, 
even if one employs the excellent ideas of Refs. [ll|, 03 • Assembling a quantum circuit from the elementary gates is 
even more difficult, particularly when the circuit is large and complicated. Unfortunately, the circuit indeed becomes 
large and complicated when one tries to apply QEC to the Hamiltonian evolution, even for the simplest case such as 
the Rabi oscillation [l(| • The largest and most complicated part of the circuit is the one that induces the encoded 
Hamiltonian evolution (such as the Rabi oscillation of a logical qubit) in a fault-tolerant manner d, ©]• Although 
a non-fault-tolerant circuit can be much simpler, such a circuit is too fragile to errors. It is therefore important to 
00 . explore new methods, which are physically more feasible and natural, for inducing the encoded Hamiltonian evolution 
' and thereby making QEC applicable. 

In this paper, we propose such a new method, choosing the Rabi oscillation as the Hamiltonian evolution to be 
protected. The method utilizes effective interactions that arise naturally among physical qubits. We show that the 
values of the parameters in the interactions are to a great extent arbitrary. Furthermore, switching of the interactions 
is unnecessary. Therefore, a system of a logical qubit with such interactions can be prepared easily by placing several 
two-level systems close to each other. Once such a system is prepared, it is driven spontaneously and flawlessly by 
the Schrodinger equation. This is much easier than to drive the system by a fault-tolerant quantum circuit. On 
the other hand, we argue physically that it is highly probable that unwanted interactions should also exist in such 
a system. While some of them are shown to be irrelevant, the others invalidate straightforward application of QEC. 
As a resolution we present a protocol, which we call the error-correction sequence. One can realize the protected 
Rabi oscillation by using the natural interactions (to induce the logical Rabi oscillation) and a quantum circuit for the 
error-correction sequence. This is much easier than realizing it wholly with a quantum circuit, because a fault-tolerant 
quantum circuit for inducing the logical Rabi oscillation, which is the largest and most complicated part of the full 
circuit, is unnecessary. 
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II. NATURAL HAMILTONIAN FOR LOGICAL RABI OSCILLATION 

We employ a two- level system as a basic element, which we call a qubit or physical qubit. We represent operators 
acting on a qubit in terms of the Pauli operators X, Y, Z (i.e., <Ti, (T 2 , (73), which are not necessarily those for a physical 
spin. To apply QEC to the Rabi oscillation, 



jXt 



|0) = cos(W)|0) + ism(ut)\l), (1) 



we replace a single qubit with a logical qubit which is composed of several qubits. The basis states |0), |1) (+1 and 
— 1 eigenstate of Z, respectively) of a qubit correspond to \0l), of a logical qubit. The subspace (of the logical 
qubit) that is spanned by the latter is called the code space. For the reasons that will be described in Scc lVIIl we 



here take the Shor code 

it 

in which a logical qubit is composed of nine qubits and 

(|000) + |111))(|000) + |111))(|000) + |111)) 
|0i) - ^ , (2) 

(|000) - |iii))(|ooo) - |iii))(|ooo) - |in)) ... 
iu) = 2375 ■ (3) 

We have to induce the logical Rabi oscillation; 

e lM \$ L ) = cos(wt)|0 L ) +isin(wt)|l z ), (4) 

where Xl is a logical Pauli operator; Xl\0l) = \1l) and Xl\1l) = \0l)- Obviously, it can be induced if the 
Hamiltonian is — loXl- [Here and after, we take h — 1.] This is an interaction among three or more qubits, for any 
code that can correct all single-qubit errors (Appendix |A"|) . For the Shor code, Xl can be represented in various ways, 
e.g., as Xl = Z^ZqZq or H^ =1 Zi, where Zi acts on qubit i. In the following, we take 

X L = (5) 

Suppose that nine qubits (such as atoms, quantum dots, and so on) composing a logical qubit are placed close to each 
other as shown in Fig.[T] Then, as will be discussed in Sec. lVIl a three-qubit interaction proportional to Xl (= Z\Z±± 
would be generated as an effective interaction. [Similar three-qubit interactions were also discussed in Refs. 
IbH].] Unfortunately, however, if this interaction is strong enough unwanted two-qubit interactions proportional to 
Z\Z^ Z^Zi, ZiZ\ should also be strong, because otherwise the following unphysical conclusion would be drawn; if one 
of qubits 1, 4, 7 is removed the other two qubits would have no interactions. Furthermore, interactions between other 
pairs of qubits, such as Z\Zi, Z2Z3, • • • , would also exist in general. Therefore, a natural and simple Hamiltonian for 
the system of Fig. [T] is 

H = H D + H s , (6) 

where 

Hd = —ujZ\Z^Zi — J{k\Z\Z^ + k^Z^Zi + kjZ^Zi), (7) 
H s = 9sZ s - X Z s . (8) 

s=2,3,5,6,8,9 

Here, u>, J, fc r 's, g s 's are real parameters. Since the signs of these parameters are irrelevant to the following discussions, 
we assume without loss of generality that they are positive. Furthermore, since three-qubit interactions are generally 
weaker than two-qubit interactions (see Sec. IVI[) . we assume naturally that 

< uj < J. (9) 

Although single-qubit terms may also exist, we can forget them because, as discussed in Appendix [Bl they are 
irrelevant to the following discussions. 

Note that the operators Z s -iZ s in Hs do not change \0l) or |1_l), i.e., they are elements of the stabilizer Q of the 
Shor code. Using this fact, we will show later by explicit calculations that the values of g s 's are irrelevant. On the 
other hand, the two-qubit interactions in Hl> are not elements of the stabilizer, and hence drive the state out of the 
code space. Nevertheless, we will show in Sec. |V] that the values of J and fc r 's are fairly arbitrary as long as ui J. 
The value of u) is also unimportant because changing oj is just equivalent to changing the time scale. Therefore, the 
values (including signs) of all the parameters in H (hence the distances between the qubits) are to a great extent 
arbitrary. This makes our scheme robust to fabrication errors. Once the system is thus fabricated, the law of the 
Nature drives it flawlessly if noise is absent. 



FIG. 1: An example of the configuration of qubits that have the proposed Hamiltonian. The distances between the qubits are 
to a large extent arbitrary. We label the qubits inside and outside the central triangle by r (= 1, 4, 7) and s (= 2, 3, 5, 6, 8, 9), 
respectively. 



III. DIFFICULTIES AND RESOLUTIONS 

We now discuss effects of noise. There are two difficulties in applying QEC straightforwardly to the system driven 
by H. We now explain them and resolutions. For simplicity, we explain the case where k\ = k± = kj = 1. More 
general cases will be discussed in Sec. |V] 

We study the first difficulty by investigating the time evolution in the absence of noise, i.e., we calculate \ip{t)) = 
e~ \ip1) , where is a vector in the code space. We note that all terms in H commute with each other, and that 
Hs does not change because all terms in Hs are elements of the stabilizer. Using these facts and the relations 
ZiZn = Z^Xl, Z±Z-j = ZiXl, Z-jZ\ = Z^Xli we find 

|^(t)) ={cos 3 (Ji) - ism 3 (Jt)} e luJtXL \iP° L ) 

+ 1 E e lJt sin(2Jt)Z r X L e^ tx -\^° L ). (10) 

r=l,4,7 

When sin(2 Jt) ^ 0, this state is out of the code space because of the last term. Therefore, we cannot perform QEC 
for phase errors at an arbitrary time, because the syndrome measurement Q to identify the errors misidentifies the 
last term as a wrong term generated by a phase-flip noise; if QEC for phase errors were performed with some intervals 
/i the time evolution would be affected as shown in Fig. ® even when noise is absent. 
However, if we focus on the discrete times 

t m = mr (m = 0,1,2,- ••), (11) 

then \ip(t m )) is in the code space, where 

t = n/2J. (12) 

Therefore, we can perform QEC at t = t m , for both phase and bit-flip errors. Furthermore, since 

m m )) = e l ^ XL K) = [cos{ut m ) + i sm(cjt m )X L ] (13) 

apart from an irrelevant phase factor, the logical Rabi oscillation is realized at these discrete times, which we call 
the discrete logical Rabi oscillation. Since to /J <C 1, the intervals r of the discrete times are much shorter than the 
period 2tt/oj of the Rabi oscillation. Hence, the discrete logical Rabi oscillation {|^(£ m ))} m =o,i,2— is quasi continuous 
as shown by the dots in Fig. [3J 

To discuss the second difficulty, let us study the time evolution in the presence of noise. Suppose, e.g., that the 
system has evolved freely from noise for t < if, where t m -i < t' < t m , until a bit-flip noise Xi acts on qubit 1 at t'. 
Then the state at t m is evaluated as 

e -iH(t m -t') Xie -iHt'^0 L) 

x(;Z 4 Z 7 re M ^ 2t '- t '"^|^). (14) 
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FIG. 2: Probability of rinding |0l) plotted against time, for the logical Rabi oscillation (chain line), the Hamiltonian evolution 
by H (dashed line), that affected by QEC for phase errors (solid line), which is performed repeatedly with some intervals p. 




Time 



FIG. 3: Probability of finding |0l) plotted against time, for the logical Rabi oscillation (chain line) and the Hamiltonian 
evolution by H (solid line). For clarity, we take u/J (<C 1) not so small; ui/J = 0.1. The dots represent the discrete logical 
Rabi oscillation. 



The terms proportional to g s 's are irrelevant because they contribute only to an overall phase factor. Therefore, g s 's 
may take arbitrary values. The problem is that the above state is different from the correctable state Xi\tp(t m )), not 
only in the term generated by Z1Z4 + ZjZ\ and {iZiZ-j) m but also in the wrong phase of the oscillation ui(2t' — t m ) - 
That is, extra errors occur because the bit-flip error in qubit 1 (or 4 or 7) is 'propagated' by H to other qubits [24| . 
As a result, QEC at t m cannot recover the correct state. 

To overcome this difficulty, we note that the syndrome measurement for bit-flip errors (unlike that for phase errors) 
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does not misidentify the state of Eq. (flQ)) as a wrong state. Hence, one can successfully perform QEC for bit-flip 
errors frequently (i.e., with intervals v which are much shorter than t) in the interval between t TO _i and t m for all m. 
As will be confirmed in the next section, this reduces the probability of errors small enough. 

Our prescription is summarized as follows: Perform QEC for both phase and bit-flip errors at all t m 's (i.e., with 
intervals r), and perform QEC for bit-flip errors repeatedly with intervals v r). The latter intervals v are not 
required to be regular. We call this protocol the error- correction sequence. 



IV. EFFECTS OF THE ERROR-CORRECTION SEQUENCE 



To see how well the error-correction sequence protects the discrete logical Rabi oscillation against noise, let us 
calculate the time evolution for to < t < t\, i.e., for < t < t, quantitatively. 

We divide the interval (0, t] into N subintervals; (0, At], (At, 2At], (2At, 3 At], • • • , where At = t/N. We model noise 
by the product of depolarizing channels [|| nf =1 £^, where £^ t acts on qubit i at the end of every subinterval as 

£®\p] = (1 - eAt)p + £ 4±it a^pa^. (15) 

Here, p denotes an input state, and e is a small positive parameter representing the strength of the interaction with 
the environment. The initial state at t — is denoted by p® L , which is assumed to be in the code space. We study its 
time evolution up to the first orders in er and lot, assuming that 

er < 1 and lot < 1, (16) 

where the latter comes from condition |9|). 

If noise and QEC were absent, p\ would evolve into 

p H (t) = e- im p° L e im = e- iHDt p° L e iHDt . (17) 

When noise is present but QEC is not performed, on the other hand, nf =1 £^] acts at the end of every subinterval. 
When N = 2, for example, p° L evolves into 



9 

'At 

t=l 



9 

e l].£ A iLPff(At)je 

i=l 

2 9 3 

e 



(18) 



= (1 - 9et)p H (t) + Ate~^^ ^ ^ *®p H (jAt)a^e iH ^ At l (19) 

j — l i—1 a—1 



By taking iV — > oo, we obtain the state at t without QEC as 

p(t,p° L )~(l-9et)p H (t) 



f alt' e -«*('-0 V V aWpH^)^^^^ - (20) 
Jo ■-, — , 



We calculate how this state is corrected by the error-correction sequence, in which bit-flip errors are corrected with 
intervals v and both bit-flip and phase errors are corrected at t = r. Although the intervals v are not required to be 
regular, and 

n = Tjv (21) 

is not required to take an integral value, we here assume for simplicity that v is regular and n is an integer. We label 
qubits in and outside the central triangle of Fig. [T] by r, r' (= 1, 4, 7) and s (= 2, 3, 5, 6, 8, 9), respectively. 

At t = QEC for bit-flip errors is performed. The pre- measurement state of the syndrome measurement is p{v, p° L ). 
The post-measurement state p'(y) depends on the outcome of the syndrome measurement. For example, when the 
bit-flip error in qubit s is detected (which happens with probability 2e^/3), 

p\v) = ^X s p H {v)X s + \y sPh {v)Y s . (22) 
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By the recovery operation, p'(v) is changed into 

p"{v) = X s p'{v)X s = \ Pa (v) + l -Z aPH {v)Z s , (23) 

which is a mixture of the correct state ph(v) and Z s ph{v)Z s , the state with a phase error in qubit s. At this stage, 
QEC for phase error is not performed because pa{ v ) is out of the code space. 

At t = 2v, QEC for bit-flip errors is performed again. The pre-measurement state is p{v, p" (v)), where p"{v) 
corresponds to one of possible outcomes of the previous syndrome measurement at t = v. We can calculate p'{2i>) and 
p"{2v) in the same way as we have calculated p'{v) and p"{v). By repeating the arguments n times, we obtain the 
probabilities of bit-flip errors during < t < t and the corresponding states p"(r) that are obtained at t = nv = t 
by correcting the bit-flip errors. To the first orders in er and lot, they are given by 



error 


probability 


corrected state p"(r) 


none 


1 - 6er (1 - 


3ct)ph(t) +f Zip H (T)Zi, 


X s 


2er/3 


\ph{t) + \Z s pb{t)Z s , 


X r 


2er/3 


\ 9 { :\t) + \Z rP { :\r)Z r , 




(t) ee f e 2UZ ^ x 
Jo 


^'p H ( T -2t')e- 2iJzM — . 



(24) 



where 

r ., „ ,dt' 

(25) 

Finally at t = r, phase errors in p"(r) are detected and corrected. We denote the state after this QEC by p"'(t). 
Since p"(r) depends on which qubit has suffered from a bit-flip error for < t < r, so does p'"{r). If a bit-flip error 
has occurred in no qubit or in qubit s, p'"(j) agrees with the correct state ph{t). If, on the other hand, a bit-flip 
error has occurred in qubit r (with probability 2er/3, see above), the conditional probability of each outcome of the 
syndrome measurement for phase errors and the corresponding p'"{r) are given by [25j ] 



error 


probability 


corrected state p"'(r) 


none or Z r 
Z r > (r' ± r) 


! + 

8 — 8 SmC lT 


KPh{t) + a^X L p H (T)X L 


+ , — 

dn ~~T~ Q"a 

±[ph(t)+X lPh (t)X l }. 



(26) 



Here, sine a; = (sina;)/x, = j§ + j^sinc^ ± |sinc^, and terms of 0(er) and 0(lut) have been dropped because 
the probability that a bit-flip error has occurred is already of 0(er). By averaging p"'(r) over all possible branches, 
we obtain the average state p c (t) under the error-correction sequence as 



Pc(t) = Ph{t) - er 



1 • 277 
1 — sine — 

n 



[ph(t) - X lPh {t)X l ] . (27) 



Therefore, p c {r) approaches the correct state Ph{ t ) with increasing n. This can be seen more clearly from their trace 
distance 0, which is calculated for n> 1 as 

itr \ Pc (t) - Ph {t)\ ~ 2n 2 eTL yz (T)/3n 2 . (28) 

Here, L yz (r) denotes the length of the projection onto the y-z plane of the Bloch vector of ph{t) in the code space. 
Hence, by taking 



n > (l/v / er)min{l,e/w}, (29) 

we can reduce the distance to about 6L yz (r) max{(er) 2 , {lot) 2 }. Since L vz {t) = 0(1), this is of the same order as 
the largest term that has been dropped in the above calculations. That is, we have successfully recovered the correct 
state at t = t (= ti), i.e., p c (ti) = Ph^i) + O (max{(er) 2 , (wt) 2 }). 

In a similar manner, we can evaluate p c (t m ) by taking p c (t m -i) as the initial state, and find that 

Pc(*m) = PH{Ua) + O (max{(er) 2 , {lot) 2 }) (30) 

for all m. Therefore, the discrete logical Rabi oscillation is protected, with only O (max{(er) 2 , [lot) 2 }) probability of 
failure, if we take n as Eq. (|2"9")l . For example, we should take n > 10 2 when er = lot — 10 -4 . 

Figure H] demonstrates how the error-correction sequence corrects errors, i.e., how the solid line approaches the 
dashed line. 



2v 



4v 



8v 



10v 



Time 



FIG. 4: Probability of finding |0r) plotted against time when p% = \0l)(0l\, for the logical Rabi oscillation (chain line), 
the Hamiltonian evolution by H (dashed line), that affected by noise (dotted line), and that corrected by the error-correction 
sequence (solid line). The dots represent the discrete logical Rabi oscillation. Here, u/J = 0.1, er = 7r/100, rjv — 6. Inset: 
magnification around t = v. 



It is clear from the results of Sees. IIIII and IIVI that the value of J is arbitrary as long as u> <C J- On the other hand, 
we have assumed in those sections that k% = = kf = 1. In this section, we show that the error-correction sequence 
is successful also when fc r 's take other values. 

Recall that the error-correction sequence consists of two parts; QEC for both phase and bit-flip errors at all t m 's, 
and QEC only for bit-flip errors with intervals v. The latter part is successful even when fc r 's are arbitrary real 
numbers, because in general a Hamiltonian which does not contain X^s and Y^s, such as the proposed H, cannot 
flip the bit of any physical qubit. Hence, the syndrome measurement for bit-flip errors does not misidentify the state 
evolved by such a Hamiltonian as a wrong state. 

Regarding the former part, we start with showing that Av's can be arbitrary integers. Note that QEC at t m 's works 
well provided that the state of the qubits at t m would be in the code space if noise were absent. As discussed in 
Sec. IIIII this condition is satisfied when k\ — k^ — k-i — 1, because \i/j(t m )) = e luJtmXL \ip1), which is certainly in the 
code space. When fc r 's are odd integers, we obtain the same result; 



V. ARBITRARINESS OF THE PARAMETERS IN H D 



m m )) 




r=l,4,7 L 



Z r Z r+3 \tp L ) 



(31) 



apart from irrelevant phase factors. Here, Z w = Z\, and we have used {ZiZ^){Z i Z'j){Z-iZi) = 1. When &v's are 
general integers (not necessarily odd), on the other hand, we have to add a certain procedure to the error correction 
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sequence. We explain this for the case where k\ = 1 and either one of £4, k 7 is even. In this case, we find that 

J>*™^|V>°) for even m, 

W mj; ~\z r Z r , e -^K) for odd m. ^ J 

Here, r and r' r) each is 1,4 or 7 depending on k^,k 7 . For example, when k^ is even and k 7 is odd, IV'(im)) = 
i?7.Z' 1 e I " t " l '' s: - L for odd m. Although this state is out of the code space, we note that the evolution into this state is 
not a stochastic process (such as evolution by noise) but a deterministic process induced by the known Hamiltonian 
H [2(|. Hence, we can surely change this state to e wtmXL \4>° L ) by applying Z 7 Zi just before QEC at t rn . By adding 
this procedure to the error correction sequence, we can successfully perform QEC at t m 's. Thus, the error-correction 
sequence, supplemented with this additional procedure, works well when fc r 's are arbitrary integers. 
Note that if fc r 's have a common factor K, one can redefine fc r 's and J as 

J' = KJ, k' r = k r /K. (33) 

The corresponding terms in Hp are then rewritten as 

k r Z r Z r +3 

= J' ^2 K Z rZ r +3- (34) 
r=l,4,7 r=l,4,7 

Hence, one can use J' instead of J, which means, e.g., that r' = tt/2J' is used instead of r. The error correction 
sequence has such flexibility. 

We next consider a more general case where Av's are rational numbers. Suppose, for example, that k\ = 1, &4 = 
3/2, /C7 = 5/3. Then, one can redefine fc r 's and J as J' — J/6,k' r — 6k r , and the corresponding terms in Hp arc 
rewritten as 

J k r Z r Z r+ a = J'(6Z 1 Z 4l + 9ZiZr + 10Z 7 Zi). (35) 

r=l,4,7 

Therefore, if one uses r' = n/2J' instead of r, the error correction sequence is successful. In general, if there exists a 
real number n such that nk r 's are integers and 

f = J/k > u, (36) 

then the error correction sequence is successful if one uses r' = tt/2J' instead of r. 

Finally, we consider the case where k^S are irrational numbers. We note that an irrational number can be well 
approximated by rational numbers. When k\ = tt (= 3.14159- • •), for example, it can be approximated by 22/7 
(= 3.14285- • ■), 333/106 (= 3.14150- ■ •), and so on. Let fci* be such a rational number. The difference k\ — kx* is 
negligible if J\k% — t<l. Therefore, for the time interval t that satisfies 

t <C 1/J\fa - fci*| , (37) 

this case reduces to the one where fc r 's are rational numbers. If one takes fci* such that \k\ — ki*\ is smaller, the upper 
limit of t given by Eq. (|37|) becomes longer, whereas condition (|36|) becomes harder to satisfy because the denominator 
of fci* becomes greater. 

To summarize this section, the error-correction sequence works well for fairly arbitrary values of fc r 's. Although it 
is better that one can successfully fabricate the system in such a way that fc r 's are integers, one can also accept most 
systems which have non- integral values of fc r 's (because of fabrication errors). This fact makes the preparation of the 
system easier. 



VI. DERIVATION OF THE EFFECTIVE INTERACTIONS 



The proposed Hamiltonian H consists of Ising-type interactions and three-qubit interactions among physical qubits. 
We here discuss how they are generated as effective interactions from more fundamental interactions. 

Many physical systems can be candidates for physical qubits that have the proposed H. As an example, we here 
consider quantum dots in a semiconductor [l5|, [l6[ . 

To be more concrete, we assume that the spin of an electron in a dot is polarized by a high external magnetic 
field, so that we can forget about the spin degrees of freedom. We also assume that the potential barrier between 
the dots is high and thick so that electron tunneling between the dots is negligible. This and (possibly) the Coulomb 
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interaction, by which states with two electrons in a single dot have much higher energies than states with a single 
electron, exclude double occupancy of a dot. For single-electron states of a dot, we assume that only the ground and 
the first excited states, denoted by |0) and |1), are relevant because higher states have much higher energies and/or 
the transition matrix elements to them are small. As a result, we can treat each dot as a system with two quantum 
levels, |0) and |1), i.e., as a qubit. For the reasons that will be explained below, we also assume that all dots in a 
logical qubit are asymmetric and different (in size and/or shape) so that accidental degeneracy is lifted. 

The effective Hamiltonian H D g for a set of such qubits (dots) is the sum of single-qubit terms and the effective 
interactions. The effective interactions are derived from more elementary interactions V, W, • • • , which are effective 
interactions among conduction electrons in homogeneous bulk semiconductors. On the other hand, V, W, • • • are 
derived from even more elementary interactions, such as the Coulomb interactions between electrons in vacuum. 
Since two- and three-body interactions have been studies in many physical systems (see, e.g., Refs. [l7], EH), we here 
consider a two-body interaction V and a three-body interaction W. Generally, the latter is much weaker than the 
former [27| . Since four- or more-body interactions are even weaker, we neglect them. 

We can represent H e s as a polynomial of the Pauli operators. In general, it would have terms that include 
Xi = \0)i + i(0| and Yi = — i|0); j(0|, where the subscript i (= 1,2, • • • ) labels the qubits. Such terms 

are non-diagonal terms that are proportional to In) (w| (m ^ n), where \n) and |m) are product states of |0)i's and 
|l)i's (such as Yli |0)i)- As discussed in Refs. [19| and [2(| and in Appendix [Cl contributions from the non-diagonal 
terms to the time evolution are negligible if 



where AE nm is the difference in energy of single qubit terms between |n) and |m). [A more precise expression of this 
condition is given in Appendix [Cl where (n\H'\m) corresponds to (n\(V + W)\m).] 

In typical situations, V and W are significant only between adjacent dots (such as dots 1,4,7, dots 1,2, dots 2,3, 
and so on, of Fig. []} because V and W generally decrease as the distance is increased. In such a case, one can 
make |A_E„ m | larger than |(n|(l/ + lF)|m)| by making the sizes and/or shapes of adjacent dots different. One can 
also make |A£Vj, m | larger by modulating spatially the magnitude of the external magnetic field. If condition (|38|) is 
satisfied by these methods, one can drop non-diagonal terms, and hence H c g reduces to H, which consists only of 
Zi = |0)ii(0| — |l)ii(l|'s, when considering the time evolution. 

On the conditions and assumptions mentioned above, H can be derived simply by taking the diagonal matrix 
elements, between |n)'s, of the effective Hamiltonian for conduction electrons, 



where Hq denotes the non-interacting part, which includes the confining potential of the dots. We here present 
explicit results for the three qubits in the central triangle of Fig. [1] Interactions between the other qubits can be 
derived more easily in a similar manner. 

Since the potential barrier is high, the wavefunctions ipr( r ) an d '4>r( r ) 01 |0) r and |l) r , respectively, are well localized 
within each dot. As a result, overlap of the wavefunctions of different dots is negligibly small, i.e., ip^{r)ijj^,{r) ~ 
for r r' and for all a, b (= 0, 1). Using this fact, we find that the effective Hamiltonian is given by 



(n\(V + W)\m) 
AE 



1 for every n, m n) 1 



(38) 



#o' + V + W, 



(39) 




(40) 



r=l,4,7 



r>r' 



where, for a, 6, c = 0, 1, 



& = c i° - 2 E(- 1 ) a (^ + v <**) - 1 E (^y^c 



(41) 




Ji4 = - jEt- 1 ) ^- - g Et- 1 )^^, 



(42) 




(43) 



a,b,c 



and similarly for the other £ r 's and J rr /'s. Here, Q? is the energy difference between |l) r and |0) r , and 




(44) 



W abc 




W(r, r', r") \^{r)\ 2 \i>\(r')\ 2 \^{r")f drdr'dr", 



(45) 
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and similarly for V a *b, V, a b- In fact, one can easily verify that all the diagonal matrix elements of Eq. ([39)l , between 
|n)'s, agree with those of Eq. (|30|) . 

It is seen that the single-dot energy ( r is renormalized by the interactions V and W, and the two-qubit effective 
interactions are generated from V and W, whereas the three-qubit effective interaction is generated from W. Regarding 
the magnitudes of the effective coupling constants, uj is much smaller than J rr <'s because the former is derived only 
from the weaker interaction W . Note that uj does not vanish by accidental degeneracy because we have assumed that 
all dots in a logical qubit are asymmetric and different. 

Since we can forget about the single-qubit terms (i.e., the first term of Eq. PO)) ') as discussed in Appendix [B] 
Eq. (140]) agrees with the proposed Hjj, Eq. ||7J), where J rr ' correspond to k r J. 

VII. DISCUSSIONS AND CONCLUSIONS 

We have shown in Sees. IIIII and IIVI that two-qubit interactions in Hjj cause errors which are correctable not by 
the straightforward application of QEC but by the error-correction sequence. One might expect that such errors 
could be corrected more easily by using more elaborate codes such as the one in Ref. [21| . If such codes are used, 
however, Xl in Hp becomes an interaction among three or more qubits. Generally, if I qubits are crowded to induce 
an ^-qubit interaction corresponding to Xl, unwanted interactions among I' (< I) qubits are also induced, as we have 
discussed on H . For any code that can correct all single-qubit errors, some of such unwanted interactions are not 
elements of the stabilizer [28[ . If I' > 3 like the code of Ref. [2l[ , they cause errors which cannot be corrected even 
by the error-correction sequence. If I' < 3 like the Shor code and the Steane code [7[, they can be dealt with the 
error-correction sequence. 

We have also shown that the values of g s 's in Hs are arbitrary. Such great flexibility would not be obtained if 
we employed a non-degenerate code because its stabilizer does not include two-fold tensor products of the Pauli 
operators. For example, the Steane code is a non-degenerate code and hence it has less flexibility. For these reasons, we 
have employed in this paper the Shor code, which is a degenerate code with I = 3 (because we can take Xl = ZiZ^Zj) 
and I' = 2. 

Possibility of use of other codes is worth exploring. It is also worth exploring the possibility of replacing a circuit for 
the syndrome measurements with another natural interactions. Our preliminary study indicates that this is basically 
possible, and more detailed studies are in progress. Furthermore, it is interesting to apply the present idea to general 
time evolutions (such as general SU(2) rotations) and/or to general systems (such as systems composed of many 
logical qubits). A possible way of realizing this may be mixed use of an Hamiltonian (such as the one of this paper) 
and simple quantum circuits. This might also be applicable to quantum simulations [22j,|23j. Since these subjects are 
beyond the scope of the present paper, we leave them as subjects of future studies. 

In conclusion, we have shown that the Rabi oscillation of a logical qubit encoded by the Shor code can be induced by 
a Hamiltonian that is composed of natural short-range interactions among physical qubits (Sec.|TT|. The Hamiltonian 
replaces the most complicated part of a quantum circuit that would be necessary for inducing and protecting the 
logical Rabi oscillation. More specifically, the state driven by the proposed Hamiltonian agrees with the logical Rabi 
oscillation at discrete times t m = mr [m = 0, 1, 2, • • •), which is quasi continuous as shown in Fig. [3] We call it the 
discrete logical Rabi oscillation (Sec. IIII[) . To prepare a physical system that has the proposed Hamiltonian, one has 
simply to place two-level systems (which are used as physical qubits), such as asymmetric quantum dots (Sec. I VI)) . as 
shown in Fig. [T] The parameters of this system, such as the positions and the sizes of the dots, are to a great extent 
arbitrary because the proposed Hamiltonian has great flexibility (Sees. ITIl and fVf . This makes the fabrication of the 
system easier. Once the fabrication is finished, one can measure the coupling constants of the effective interactions, 
and the important parameters such as r can be calculated from them. To protect the discrete logical Rabi oscillation 
against noise, the ordinary QEC cannot be applied straightforwardly. However, we have shown that it can be protected 
by a new protocol, which we call the error-correction sequence (Sees. IIIII and ITV|) . In this protocol, QEC for both 
phase and bit-flip errors is performed at t TO 's, whereas QEC only for bit-flip errors is performed frequently in the 
interval between t m _i and t m for all m. The frequency of the latter is determined by the strength of noise and the 
parameters of the effective interactions (Sec. IIV[) . One can realize the protected Rabi oscillation by using the natural 
Hamiltonian (to induce the logical Rabi oscillation) and a quantum circuit for the error-correction sequence. This is 
much easier than realizing it wholly with a fault-tolerant quantum circuit. 
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APPENDIX A: X L IS AN INTERACTION AMONG THREE OR MORE QUBITS 

Let P c be the projection operator onto the code space; 

Pc = |0L)<0 L | + |l i )(l i |. (Al) 
An rt-qubit code which can correct all single-qubit errors satisfies the following condition Q; 

P c a$afp c = Xla , 3 pPc (i,j = l,--- ,n; a,(3 = 0,1,2,3). (A2) 

Here, a$ denotes the identity (a = 0) and Pauli (a = 1, 2, 3) operators acting on qubit i, and Xia,jp i s an element of 
some Hermitian matrix. 

If Xl were a Pauli operator or a two-fold tensor product of Pauli operators, Eq. (|A2|) could not be satisfied. For 
example, if Xl = X\X 2 for some code the left-hand side of Eq. (|A2j> with a„ = X\, = X 2 (i.e., i = l,j = 2, a = 
[3 = 1) reduces to 

P C X X X 2 P C = P C X L P C = |0l><1l| + |U>(0l|. (A3) 

Since this is neither vanishing nor proportional to P c , Eq. (|A2|) is not satisfied for any value of xii- This means that 
such a code cannot correct all single-qubit errors. 

Therefore, Xl is a three- or more-fold tensor product of the Pauli operators (which corresponds to an interaction 
among three or more qubits) for any code that can correct all single-qubit errors. 



APPENDIX B: IRRELEVANCE OF SINGLE-QUBIT TERMS 

When two levels of a qubit have different energies, a single-qubit term, which represents the energy difference, arises 
in its effective Hamiltonian as discussed in Sec. IVI1 All effects of such single-qubit terms can be canceled if we do 
everything in the rotating frame (20| . Although this fact seems to be known widely, we here explain it in order to 
clarify its meaning in the context of QEC. 

Let us investigate the time evolution of a state p+ by the following Hamiltonian 

9 9 

H+ = - - ]T &Zi + H = --J2 CiZi + H D + H s , (Bl) 

i=l i=l 

where Q's are real numbers. We can go to the rotating frame (an interaction picture) by Uq = exp(^ X^j=i Citify, as 
p mt = UqP+Uo. It evolves according to 

ij/ ot - PtHUo, / ot ] = [H, p ro % (B2) 

where we have used [H, Uo] — 0. Therefore, p rot undergoes the same unitary evolution as that of p of Sec. IIVI 
Furthermore, it is easy to show that the depolarizing channel in the rotating frame is also the same as the one in 
Sec. IIVI Thus, in the presence of noise, p mt evolves in the same manner as p of Sec. IIVI Therefore, the error correction 
sequence will be successful if we set the initial state p lot (0) in the code space and perform QEC in the rotating frame. 
For example, the observables for the syndrome measurement in the rotating frame are M£°* = Z1Z2, = 

X1X2X3X4X5XQ, and so on. In the laboratory frame (Schrodinger picture), they are given by M bl = U ZiZ 2 U$ = 
Z ± Z 2 and M Pl = UoX^X^XiX^ul = lf l=1 e^-p(^QZ l t)X l exp(-^QZ l t), respectively 



APPENDIX C: IRRELEVANCE OF TERMS INCLUDING X, Y 

It seems widely accepted by researchers of NMR that the non-diagonal terms, which include Xi's and/or Y^s, in 
H e g are irrelevant to the time evolution if condition (|38|) is satisfied (see, e.g., Refs. [lj| and [13|). For completeness, 
we here show that this is indeed true under reasonable assumptions. 

Let us decompose H e s as 



H cB = H + H + H\ H =-^J2 ^ 

i 



(Cl) 
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where Q is the energy difference (that is renormalized, like Eq. (|4Tj) . by interections among dots) between |1)^ and 
\0)i, H is a polynomial of Z^s only, and H' consists of the other terms (such as X\Yi, XiX 2 Z 3 , and so on) which 
include Xi's and/or Y^s. 

We denote a product state of |l)i's and |0)j's, such as FL by I")- ln terms of such product states, Hq and H 
are diagonal, whereas H' gives the off-diagonal elements. To characterize the magnitude of the latter, we define the 
parameter £„ m by 

(n\H'\m) , , 

it (n\H \m) f= 0, 



£nm = ^ Ai?„ m (C2) 
if (n|fT'|m> = 0, 

where AE nm denotes the difference of the eigenvalues of H between \n) and \m). We also define 



f=w£i&»»i a - ( C3 ) 

Consider the time evolution operator U e s(t) generated by H c r. We can write it as 

U cS (t) = cxp (-iflefft) = [/(f)e" iQ W , (C4) 

where 

£/(t) = e -*(«o+H)i = e -iH t e -iSt f ( C5 ) 

and is the Hermitian operator defined by e~ % ^' = W(t)U e s(t). It is clear that 

Q{t) = when £ = 0. (C6) 

If ^ were large then Q(t) would be significant, particularly when AE nrn = for all n,m, for which £ = oo. On the 
other hand, if £ is small enough ||Q(i)|| is expected to be small, where || • || denotes the operator norm. It is natural 
to assume that 

Assumption 1: Q(t) is continuous in £ nm 's in the neighborhood of £ = 0. (C7) 

This assumption seems reasonable from the perturbation expansion of the time evolution operator in the interaction 
picture, which corresponds to e lHat U e fi{t) = e~ lHt e~ l Q^\ 

l-itJ2(n\H\n)\n)(n\-i Y, I e lAjE ""**' dt'(n\H'\m)\n)(m\ + ■ ■ ■ 

= l-itJ2in\H\n)\n)(n\-J2 E (e lAEnmt - l) £ n m\n){m\ H , (C8) 

n n m (=£ n ) 

each term of which is continuous with respect to £ nm . 

Assumption 1, together with Eq. (|C6j) . means that for any small positive number e there exists a positive number 
£ E) t such that 

\\Q(t)\\ < e for all £ < ^, t . (C9) 

In other words, for a given time period [0,t) we can neglect Q(t), i.e., we can regard U e g(t) = U(t), if £ is small 
enough. This means that the time evolution by H e g takes place as if |n)'s (which are eigenstates of H + H) were its 
eigenstates. That is, if we expand an initial state in terms of |n)'s as ^ c n \n), 

e-^t^c^n) ~J2cne-* m ' +H)ln)t \n). (CIO) 

n n 

Note that the above argument is general in the sense that we have not assumed any specific forms for H and H'. 
For example, the argument in Appendix A of Ref. [l9j |. where specific forms have been assumed, is essentially a special 
case of the present general argument. 
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In the above argument, we have not excluded the possibility that £ e ,i increases with increasing t. This will not 
cause difficulty when one sets an upper limit of t. To be more complete, however, we here discuss dependence of £ e> t 
on t. We note that the coefficients of the third term of Eq. 1|C8]) are upper bounded as 

| ( e iAE„ m t _ y < 2 |^ m | < ^ (cu) 

for all t. This is due to the fact that t appears only through the oscillatory factor e' AE ™'. Since this is the case also 
for higher-order terms, we expect that 

Assumption 2: £ £jt has an upper bound £ e , which is independent of t. (C12) 

If this is true, then for any small positive number e and for all t 

\\Q{t)\\ < £ for all£<f £ . (C13) 

In other words, we can regard U e s{t) = U(t) even for long t if £ is small enough. 
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